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was  quasi-steady  for  a current  discharge  period  of  800  microseconds.  Gas  injecti< 
was  continuous  and  was  dona  in  two  ways  - in  one  mode  it  was  introduced  through 
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the  magnet  from  piano  wires  and  measuring  its  deflection  when  the  thruster  was 
fired.  For  argon  injection  through  the  electrodes,  the  current  was  carried  down- 
stream from  the  magnet  in  the  usual  self-field  discharge  mode  and  no  magnet 
deflection  occurred.  For  axial  injection  of  argon  a large  magnet  deflection 
occurred.  Heweveic,  ^his  was  found  to  be  due  to  current  carried  on  the  sidewalls 
which  was  stabilized  by  the  sidewalls  and  which  produced  severe  ablation  of  the 
alumina.  Whan  a mixture  of  hydrogen  and  nitrogen  was  substituted  for  argon,  no 
magnet  deflection  occurred^^JLndicating  the  absence  of  wall  current  and  presumably 
ablation.  In  this  test  a ceramic  rod  spanning  the  channel  from  one  electrode  tc 
the  other  to  stabilise  the  current  in  the  magnet  gap  (flame  holder)  had  no  affect 
It  was  concluded  that  impulse  enhancement  through  the  use  of  a bias  magnetic 
field  is  not  feasible. 
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1 . INTRODUCTION 
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1.1  RELEVANCE  TO  AIR  FORCE  NEEDS 

Low  total  impulse  requirements  for  Air  Force  satellites  are  presently 
provided  by  small  monopropellant  rocket  engines  and  resistojets.  The 
growth  in  total  impulse  capability  needed  for  the  near  future  is  being  met 
by  improving  the  useful  lifetime  of  these  engines.  However,  future  satel- 
lites requiring  much  larger  total  impulse,  either  for  attitude  control  and 
station  keeping  for  long  periods  or  for  larger  impulse  maneuvers,  such  as 
orbit  repositioning,  will  require  the  high  specific  impulse  available  with 
electric  propulsion  to  avoid  the  unacceptable  increases  in  propellant  mass 
of  the  lower  I chemical  thrusters.  Development  of  the  pulsed  ablative, 
solid  propellant  thruster  is  proceeding  for  the  relatively  low  impulse 
stationkeeping  missions.  To  perform  orbit  repositioning  with  lighter 
weight  than  the  chemical  systems  now  in  use,  and  to  provide  the  capability 
for  future  requirements  of  multiple  orbit  shift  maneuvers,  a compact, 
lightweight,  reliaole  electric  thruster  is  needed.  A pulsed  electro- 
magnetic thruster  employing  a strong  bias  magnetic  fie^  appeared  to  be 
a candidate  for  meeting  this  requirement  with  high  efficiency. 

1.2  MAGNETIC  FIELD  ENHANCEMENT 

Plasma  thrusters  use  the  electromagnetic  force  developed  by  the  plasma 
current  in  the  local  magnetic  field  to  accelerate  the  plasma  mass  away  from 
the  thruster.  In  most  cases  the  magnetic  field  is  the  resultar-t  of  the 
fields  due  to  current  flowing  in  the  circuit  and  in  the  plasma,  i.e.,  the 
"self-fields."  For  self-field  accelerators,  the  current  density  at  any 
point  in  the  plasma  is 

j - o(E  - vB)  (1) 

where  g is  conductivity,,  v is  plasma  velocity,  and  E and  B are  electric 
and  magnetic  fields.  The  directions  of  v,  E,  and  B are  mutually  orthogonal 
and  j is  parallel  to  E. 


Equation  (1)  can  be  rearranged  to  give 


E * ^ + vB 


(21 


After  multiplication  by  j this  becomes 

,2 

jE  - * (JB)v 


(3) 


which  expresses  the  Input  power  density  as  the  sum  of  the  ohmic  dissipation 
power  density  and  the  rate  of  mechanical  work  per  unit  volume  done  by  the 
jB  force  density.  If  the  only  losses  were  ohmic,  the  efficiency  would  there- 
fore be 


JB 


or 


jBv  +_j 
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1 

TT 


(4) 
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For  the  simple  case  of  a slab  model  with  no  bias  field,  with  property 
variations  In  the  acceleration  direction  only,  the  magnetic  field  gradient 
is  approximately  B/5  where  6 is  the  Interaction  distance.  Thus  the 
Maxwell  Curl  - H relation  Is 


J V 6 


where  u Is  the  permeability  of  space,  and  Equation  (4)  can  be  written 


(5) 


1 + mr 


UOJV 


(6) 


but  since  the  diffusion  speed  of  the  magnetic  field  relative  to  the 
plasma  is 


v 


d 


1 

U<j<5 


(7) 


Equation  (6)  becomes 


n 


(8) 


where  the  magnetic  Reynolds  number  R^  is  the  ratio  of  the  plasma  velocity 
to  the  diffusion  speed.  By  multiplying  the  velocity  terms  by  the  accelera- 
tion time  it  Is  seen  that  the  Reynolds  number  Is  also  equal  to  the  ratio 
of  the  acceleration  distance  to  the  field  penetration  depth. 


R 


N 


(9) 


This  ratio  Is  large  compared  to  unity  for  electromagnetic  machines  because 
the  magnetic  force  Is  applied  over  a distance  much  larger  than  the  dimen- 
sion of  the  winding.  This  is  also  true  for  a metal  plate  accelerated  impul- 
sively by  the  rising  current  in  a coil  placed  close  to  the  plate,  and  for 
an  exploding  supernova. 

For  attitude  control  and  station  keeping  of  present  satellites,  the 
acceleration  distance  is  comparable  to  the  current  sheet  thickness  so  that 
the  Reynolds  number  Is  nearly  unity  and  the  efficiency  is  limited  to  about 
50X  because  of  ohmic  losses  alone.  Also  for  steady  or  quasi-steady  plasma 
thrusters  the  current  and  magnetic  field  penetrate  the  plasma  for  the  same 
distance  over  which  acceleration  occurs  and  the  efficiency  for  such  devices 
is  also  limited  to  50X. 


This  limitation  was  recognized  a few  years  ago  during  a study  of 
Impulsive  plasma  acceleration  by  inductive  coupling  between  the  plasma  and 
circuit  currents  that  was  conducted  by  TRW  Systems  for  the  Air  Force  Office 


of  Scientific  Research.  An  attempt  was  made  at  that  time  to  Increase  the 
magnetic  field  in  the  Interaction  region  through  the  use  of  inductively 
produced  bias  magnetic  fields.  The  first  experiment  of  this  kind  was 
successful.  It  used  a pair  of  bias  field  coils  whose  axes  were  ortho- 
gonal to  the  accelerator  coil  axes  so  that  the  bias  and  accelerator  cir- 
cuits were  decoupled.  Attempts  were  then  made  to  incorporate  this  prin- 
ciple in  a practical  thruster  configuration  by  using  coaxial  bias  and 
accelerator  coils.  Shorted  turn  shield  coils  were  used  to  decouple  the 
circuits.  It  was  found  that  fields  due  to  the  current  pulses  in  the  shield 
coils  effectively  eliminated  the  bias  fields  at  the  time  the  plasma  cur- 
rent was  developed.! 

The  program  described  in  this  report  was  an  attempt  to  avoid  the 
limitation  of  the  inductive  bias  field  accelerator  by  using  a permanent 
magnet  to  produce  the  bias  field  in  an  electrode  thruster.  It  should  be 
mentioned  that  a large  scale  impulsive  thruster  offers  another  way  to  get 
around  the  low  Reynolds  number  limitation,  provided  "slug  model"  accelera- 
tion can  be  used  to  avoid  the  50%  impact  heating  loss  that  occurs  in  "snow 
plow  model"  acceleration.  This  may  be  feasible  for  future  large  scale 
satellites. 


^ Dally,  C.L.,  H.A.  Davis,  and  R.H.  Lovberg,  "Magnetic  Field  Annihilation 
of  Impulsive  Current  Sheets,"  AFOSR-TR -76-0698,  March  1976. 
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2.  THEORETICAL  BACKGROUND 

The  physical  arrangement  of  the  bias  field  thruster  is  indicated 
schematically  in  Figure  1.  The  plasma  is  accelerated  in  a V-shaped  ceramic 
channel  that  is  placed  between  the  pole  pieces  of  a permanent  magnet.  Cur- 
rent flows  through  the  plasma  between  the  upper  and  lower  electrodes.  The 
thruster  operates  in  the  quasi-steady  mode  with  gas  entering  the  channel 
through  a narrow  sonic  slit  and  accelerating  to  a pressure  less  than  1 torr 
in  the  current  carrying  region.  In  addition  to  accelerating  the  gas,  the 
divergent  channel  also  serves  to  decrease  the  magnetic  field  In  the  acceler- 
ation direction.  The  reason  for  this  can  be  seen  by  considering  first  the 
efficiency  obtainable  in  a constant  area  channel. 


GAS  INLET 


WALLS 


Figure  1.  Schematic  Arrangement  of  Bias 
Field  Thruster 


2.1  CONSTANT  BIAS  FIELD 

It  Is  Interesting  that  the  application  of  a bias  field  that  Is  both 
spatially  and  temporally  constant  does  not  alter  the  ohmic  efficiency.  At 
any  location  In  the  accelerator 


0V$-  » 

OI-  _dv..JB  (,0) 

dZ  T 

where  p Is  mass  density,  Z Is  the  acceleration  direction  and  m Is  mass  flow 
rate  per  unit  area.  For  steady  flow,  and  for  a bias  field  that  Is  large 
compared  to  the  self-field,  both  m and  B are  constant.  When  Equation  (1) 
is  differentiated 


- oB 


dv 

HZ 


or  with  Equation  (10) 
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which  Integrates  to 
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(11) 
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This  relation  for  the  current  density  can  be  combined  with  Equation  (10) 
to  give 


4-* 


which  Integrates  to 


.A 


#v- 


-4-a 


(12) 


From  Equation  (11)  the  ohmic  dissipation  power  density  Is 

>2 


p -4  ■¥ 


-2a  B 


(13) 


and  by  combining  Equations  (11)  and  (12)  the  kinetic  power  density  is 


< 2 ( -oB?Z  2 q B2z\ 
■ JvB  \e  * - e ifi  / 


(14) 


By  Integrating  Equations  (13)  and  (14)  It  Is  found  that 
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The  variations  with  Z of  the  dlsslpattve  and  kinetic  power  densities 
are  shown  In  Figure  2.  The  former  Is  a maximum  at  Z * 0 and  drops  one 
e-fold  by  the  time  the  kinetic  power  reaches  Its  maximum.  It  Is  seen  that 
not  only  B,  but  a and  m have  no  effect  on  the  50%  limit  on  efficiency.  The 
e-fold  distance  is  aB2/m  and  Increasing  the  bias  field  with  all  other  vari- 
ables held  constant  simply  reduces  the  scale  of  the  acceleration  distance 
without  affecting  the  efficiency. 
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Figure  2.  Relative  Dissipation  and  Kinetic  Power 
Densities  (Constant  E and  B) 

2.2  DECREASING  BIAS  FIELD 

The  ohmic  efficiency  limit  of  50%  Is  not  Inherent  In  the  bias  field 
concept  but  Is  a consequence  of  the  bias  field  being  held  constant.  It  can 
be  removed  by  using  a field  that  decreases  In  the  acceleration  direction. 

The  reason  for  this  can  be  seen  Intuitively  by  referring  to  Equation  (1). 

The  current  density  (and  hence  the  acceleration)  drops  to  zero  whenever 
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v * E/B , l.e.,  when  the  back  emf  (In  electric  motor  language)  t$  equal 
to  the  applied  electric  field.  This  velocity  (the  plasma  drift  speed)  is 
low  when  B is  large.  Thus  by  applying  a large  bias  field  near  Z * 0 the 
50%  limit  on  ohmic  efficiency  is  reached  when  v Is  small,  i.e.,  when  both 
the  kinetic  and  ohmic  energies  are  low.  The  B field  is  then  decreased,  at 
constant  E,  so  that  E/B  Increases  thereby  increasing  the  kinetic  energy. 

This  behavior  can  be  examined  In  more  detail  by  assuming  that  the 
bias  field  varies  Inversely  as  the  local  pole  face  separation.  A linearly 
divergent  channel  is  used,  as  illustrated  in  Figure  3. 


Figure  3.  Accelerator  Channel  Geometry  for  Decreasing 
Bias  Field 

The  distance  from  the  centerline  to  the  channel  wall  is  y (at  any  Z) 
and  the  distance  to  the  magnet  pole  face  is  yg.  If  Bq  represents  the  field 
strength  at  Z a 0,  then  the  value  at  any  Z is 


but  since 

y8  1 yoB  * Z “ 
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Another  consequence  of  the  divergent  channel  geometry  is  that  m Is  no 
longer  constant.  If  ft  represents  the  total  mass  flow  rate  and  h is  the 
channel  height  then 
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The  circuit  current  is  the  integral  of  the  current  density  over  the 
electrode  length  L, 


L L 

I » 2 / jydZ  » 2y0  /(I  +i-Tana)jdZ 
0 0 0 


(19) 


and  the  local  plasma  density  is  given  by 


. , M . M 

**  Vv  (1  .-Z-Tana) 
yo 


(20) 


A numerical  procedure  can  be  used  to  calculate  accerator  efficiency 
by  using  Equation  (18)  to  calculate  j at  an  axial  position  Z where  E,  BQ 
and  v are  known.  Equation  (17)  then  gives  dv/dZ  and  the  velocity  at  an 
adjacent  downstream  station  becomes 


v + Av  » v + (21 ) 

2.2.1  Dimensionless  Formulation 

If  the  current  density  from  Equation  (1)  Is  substituted  into  Equa- 
tion (10),  the  acceleration  relation  becomes 


dv 

dZ 


^-(E-vB) 


or 


dv 

dZ 


oB" 


p<e7b> 


■d 


(22) 


where  E/B  is  the  local  drift  speed.  By  making  the  Initial  value  of  B large 
(where  the  magnet  gap  Is  a minimum)  the  plasma  speed  approaches  the  drift 
speed  in  less  than  a millimeter.  The  remainder  of  the  acceleration  occurs 
with  a small  value  for  the  difference  term  in  Equation  (22)  multiplied  by 
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a large  coefficient.  The  magnitude  of  this  coeffclent  Increases  with  the 
conductivity,  with  the  square  of  the  local  bias  field  and  with  decreasing 
local  density.  This  Illustrates  how  a large  bias  field  can  compensate  for 
a low  conductivity. 

This  point  can  be  seen  In  a different  way  by  completing  the  dimension- 
less formulation  of  Equation  (22).  Thus 


oL  B M v x 
~pv~  (1  - OT 


where  L Is  accelerator  length,  or  since 


q L B 
pv 


It  follows  that 


sAi. 


The  Aflven  Mach  number  MA  Is  the  ratio  of  the  plasma  velocity  to  the  propagation 
speed  of  magnetic  waves  In  the  acceleration  direction.  By  making  this  number 
small,  through  the  use  of  a strong  bias  field,  the  effective  Reynolds  number 
Is  Increased. 

A more  useful  dimensionless  formulation  for  calculation  Is  to  refer 
the  velocity  to  the  exit  value  of  the  drift  speed  rather  than  to  the  local 
value.  Equation  (23)  then  becomes 
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The  values  for  RN  and  MA  are  found  after  a calculation  Is  finished  and  <*E 
has  been  determined.  Thus  the  exit  values  of  and  are  given  by 
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From  Equation  (1)  the  dimensionless  current  density  Is 


The  circuit  current  is 
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The  Input  power  Is 
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and  the  ohmic  loss  is 
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The  ohmic  efficiency  is  therefore 
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From  Equation  (29)  a qualitative  estimate  of  the  efficiency  is 
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or,  from  Equation  (24) 


dn  RN 
MA 


0 - n ) 


After  integration  this  becomes 
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A plot  of  this  relation  is  shown  in  Figure  4.  The  efficiency  Is  over- 
estimated but  the  importance  of  the  Reynolds  number  parameter  is  correctly 
il  lustrated. 


Figure  4.  Effect  of  RN/M^  on  Ohmic  Efficiency 

2,  TT  “ 10.  « » 20°) 


r -•  1 


2.2.2  Effect  of  Accelerator  Parameters  on  Ohmic  Efficiency 

Figure  4 shows  the  effect  of  the  parameter  RN/M^  on  ohmic  efficiency 
for  the  geometric  parameters  yoB/yQ  * 2,  L/yQ  » 10  and  a - 20°.  At  a value 
of  1.0  the  efficiency  Is  40%,  which  agrees  well  with  the  self-field  acceler- 
ator estimate.  At  a value  of  10  the  efficiency  has  Increased  to  84%. 

The  corresponding  accelerator  parameters  are  plotted  In  Figure  5. 


The  effects  on  efficiency  of  the  magnet  to  channel  gap  ratio,  the 
channel  length  to  gap  ratio  and  the  channel  wall  angle  are  shown  In  Fig- 
ures 6 through  8,  respectively.  From  Figure  6 It  is  seen  that  the  gap  ratio 
should  be  as  small  as  possible  for  maximum  efficiency.  This  means  that  the 
Insulating  walls  separating  the  magnet  pole  faces  from  the  plasma  should  be 
as  thin  as  possible.  Figure  7 shows  that  the  channel  length  should  be 
about  10  times  the  gap  dimension. 

The  optimum  length  would  require  a tradeoff  between  efficiency  and 
magnet  weight,  which  Increases  with  channel  length. 
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Figure  5.  Dimensionless  Accelerator  Parameters  vs 
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3.  DESCRIPTION  OF  ACCELERATOR 


3.1  CHANNEL 


Several  accelerator  channels  were  tested  during  the  program  In  an 
attempt  to  decrease  the  amount  of  Insulator  oblation  that  was  added  to  the 
plasma  during  the  discharge.  The  Initial  channel  Is  shown  In  Figure  9. 

It  was  mounted  to  a luclte  plate  on  one  of  the  ports  of  a 4 x 8 foot 
vacuum  chamber  as  shown  In  Figure  10. 

The  electrodes  of  the  channel  shown  In  Figure  9 were  made  of  l/4-1nch 
thick  brass  while  1/16-1 nch  thick  glass  plates  were  used  for  the  side  walls. 
Gas  was  admitted  to  the  acceleration  region  through  a sonic  slit  about 
1/4  mm  wide.  The  electrode  separation  was  2.5  cm.  The  electrode  length 
In  the  flow  direction  was  3.2  cm  for  the  first  device  tested.  This  was 
later  shortened  to  1.3  cm  by  covering  the  downstream  portions  with  glass 
plates.  The  last  channel  tested  Is  shown  In  Figure  11.  It  had  tungsten 
electrodes  and  all  four  of  the  channel  surfaces  were  made  of  high  purity 
alumina.  The  edges  where  the  surfaces  joined  were  recessed  such  that  the 
RTV  used  for  the  vacuum  seal,  seen  In  Figure  11,  was  not  exposed  to  the 


Figure  10.  Accelerator  Channel  Mounted  on  Vacuum  Chamber 


plasma.  For  this  model  the  only  materials  subject  to  ablation  were  alumina 
and  tungsten. 

3.2  MAGNET 

Recently  developed  Brown-Bovarl  samarium-cobalt  permanent  magnets  were 
used  to  drive  the  bias  field  magnetic  circuit.  A demagnetization  curve  for 
this  material  Is  shown  %1n  Figure  12.  As  shown  in  Figure  12,  eight  magnet 
blocks  were  used  in  the  structure  (four  in  parallel  and  two  in  series). 

Each  of  the  blocks  was  1.25  inches  long  and  had  a 1-inch  square  cross 
section. 

The  original  design  used  four  magnetic  blocks  with  a minumum  gap  spac- 
ing of  0.19  cm.  The  field  measurements  plotted  in  Figure  14  show  that  high 
fields  were  obtained  for  the  first  few  millimeters  with  low  values  near  the 
exit.  The  gap  was  opened  to  0.51  cm  to  provide  enough  room  for  the  glass 
side  walls  so  that  the  apex  of  the  plasma  volume  would  be  near  the  leading 
edge  of  the  pole  pieces.  Field  strength  plots  are  shown  in  Figure  14  for 
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Figure  11.  Photograph  of  Ceramic  Channel  Showing  Gas 
Manifold  and  Current  Return  Lead 

four,  six,  and  eight  magnet  blocks.  The  latter  configuration  was  used  for 
most  of  the  experiments.  As  shown  in  the  figure,  the  field  along  the  axis 
is  prescribed  reasonably  well  represented  by  the  relation 

r * 0-954 

By  . . 7z  Tan  IZ.V* 

1 ‘ 0.51  - 

The  lateral  variation  of  the  magnetic  field  across  the  gap  was  measured 
at  several  positions  to  determine  the  uniformity  of  the  field  at  a given 
axial  station.  The  results,  plotted  in  Figure  15,  show  that  the  variation 
was  reasonably  small. 

For  the  last  series  of  tests  conducted  during  the  program,  those  con- 
ducted with  the  ceramic  accelerator  channel  shown  in  Figure  11,  the  eight 
magnet  blocks  were  remounted  in  the  pole  pieces  shown  in  Figure  16.  For 
this  magnet,  the  narrow  end  gap  was  increased  to  0.35  inch  and  the  length 
of  the  pole  face  (in  the  acceleration  direction)  was  reduced  to  0.7  inch. 


Figure  13.  Bias  Field  Magnet 


This  moved  the  magnet  forward  relative  to  the  electrodes  so  as  to  increase 
the  bias  field  in  the  plasma  current  region. 


3.3  CURRENT  SOURCE 


Current  for  the  discharge  was  supplied  by  the  six  element  pulse  line 
shown  in  Figure  17.  Each  section  of  the  line  had  780  pF  capacitance  and 
6 pF  inductance.  The  total  capacitance  of  the  electrolytic  capacitors 
was  4800  pF  and  the  pulse  duration  of  the  line  was 


Figure  15.  Variation  of  Magnetic  Field  Across  Magnet  Gap  for 
Several  Axial  Stations  (eight  magnets,  6 = 0.51  cm) 


Figure  16.  Large  Gap  Magnet  Used  with  Ceramic  Channel 
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4.  EXPERIMENTAL  METHODS 


4.1  EFFICIENCY  MEASUREMENT 

The  experimental  efficiency  was  found  from  the  relation 
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where AI  is  the  impulse  per  shot,  M is  propellant  ejected  per  shot  and  J 
is  the  product  of  voltage  and  current  at  the  electrode  terminals  integrated 
over  the  discharge  period.  The  I$p  was  found  from  the  relation 
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The  mass  flow  rate  Into  the  thruster  was  calculated  using  the  measured 
fillup  rate  of  the  vacuum  chamber  at  constant  supply  pressure. 


The  impulse  bltAl  was  determined  by  suspending  the  magnet  structure 
from  a pair  of  1.4  meter  long  piano  wires  and  measuring  its  deflection  when 
the  thruster  was  fired;  the  deflection  was  of  the  order  of  2 mm  and  could 
be  read  with  a repeatibility  of  ±5*  by  using  a microscope. 


The  impulse  measured  in  this  way  included  only  the  electromagnetic 
force  arising  from  the  interaction  of  the  plasma  current  with  the  bias 
field;  it  did  not  include  the  impulse  due  to  the  self-field  acceleration. 
Thus  zero  magnet  deflection  would  indicate  no  gain  in  efficiency  relative 
to  a self-field  accelerator,  rather  than  zero  thruster  efficiency. 


An  independent  measurement  of  total  impulse  was  made  by  collecting  the 
accelerated  plasma  In  a 12-Inch  diameter  by  16-inch  long  cylinder  of  alu- 
minum foil  suspended  from  threads  Inside  the  vacuum  chamber.  Four  threads 
were  used,  attached  to  a lightweight  frame  to  which  the  cylinder  was 
mounted.  This  method  of  suspension  minimized  the  oscillations  that  were 
Induced  by  the  high  velocity  gas  flow  through  the  accelerator  before  the 
capacitor  bank  was  discharged.  The  data  reading  technique  was  first  to 
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observe  the  position  of  a pointer  attached  to  the  cylinder  relative  to  a 
meter  stick  in  the  vacuum  chamber  with  no  gas  flow,  then  to  estimate  the 
maximum  deflection  in  the  gas  stream  before  and  after  the  discharge.  From 
the  known  geometry  of  the  suspension  the  increase  in  potential  energy  of 
the  pendulum  was  calculated. 

The  input  energy  in  the  denominator  of  Equation  (33)  was  measured  by 
integrating  the  product  of  the  circuit  current  and  the  voltage  across  the 
electrodes  for  each  test  condition.  A typical  oscillograph  record  is  shown 
in  Figure  18.  A pair  of  voltage  probes  and  a differential  amplifier  were 
used  for  the  voltage  data  while  the  voltage  drop  across  a low  impedance 
series  resistor  (0.167nn.')  was  used  to  measure  current.  The  accelerator 
impedance  for  this  case  is  seen  to  be  about  50  nn.'. 
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Figure  18.  Current  (Upper  Trace)  and  Terminal  Voltage  (Lower  Trace) 
for  a Typical  Discharge  (30  V/cm,  1200  A/cm) 

4.2  DIAGNOSTIC  MEASUREMENTS 

Measurements  of  current  and  magnetic  fields  in  the  plasma  were  made 
with  miniature  probes.  Because  of  the  small  sire  of  the  thruster  it  was 
not  feasible  to  obtain  meaningful  spatial  distribution  of  these  properties. 
The  measurements  however  were  helpful  in  locating  the  current,  i.e.,  in 
determining  the  principal  features  of  the  discharge. 
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5.  RESULTS 


As  can  be  seen  from  Equation  (32)  the  ratio  ofnto  Isp  Is  Independent 
of  the  mass  accelerated  and  depends  only  on  measurements  of  Al  and  J. 
Measurements  with  the  first  accelerator  tested  (Figure  9)  gave  I$p  values 
calculated  from  Equation  (33)  that  led  to  efficiencies  from  Equation  (32) 
well  over  100%.  There  appeared  to  be  two  possible  explanations.  One  was 
that  the  current  ran  partly  along  the  glass  sidewall  and  the  other  was  that 
additional  mass  was  added  by  ablation.  If  current  were  running  on  the  side- 
wall  it  would  produce  a large  force  on  the  magnet  without  exerting  a com- 
parable force  on  the  plasma  because  of  wall  friction.  To  examine  this 
possibility  the  lateral  distribution  of  the  axial  magnetic  field  was 
measured  midway  between  the  electrodes.  If  this  variation  were  monotonic 
it  would  indicate  current  flow  on  the  sidewalls,  while  if  maxima  and 
minima  were  found  it  would  show  that  the  current  was  confined  between  those 
peaks. 

As  shown  in  Figure  19,  the  axial  magnetic  field  distribution  at 
2.9  cm  from  the  gas  injection  slit  showed  well  defined  maxima  and  minima. 
These  features  were  not  observed  at  1.6  cm  and  it  is  possible  that  current 
was  flowing  on  the  glass  sidewalls  at  that  location.  Figure  20  shows  that 
a sharply  peaked  maximum  in  the  axial  variation  of  the  same  magnetic  field 
component  occurred  close  to  the  2.9  cm  position.  Since  this  corresponded 
to  the  location  of  the  epoxy-fiberglass  plate  to  which  the  channel  was  at- 
tached (Figure  9)  it  was  felt  that  ablation  was  the  cause  of  the  spuriously 
high  values  of  the  measured  impulse. 

Inspection  showed  extensive  charring  of  the  plate  had  occured.  To 
check  this  possibility  further,  the  plasma  momentum  was  measured  with  the 
ballistic  pendulum  suspended  inside  the  vacuum  chamber.  The  result,  shown 
in  Figure  21,  indicates  good  agreement  between  the  magnet  impulse  and  the 
plasma  impulse.  This  seemed  to  show  that  there  could  not  have  been 
appreciable  current  flowing  on  the  glass  sidewalls. 

Four  glass  plates  were  attached  to  the  inner  surfaces  of  the 
accelerator  that  covered  the  plate  and  electrodes  except  for  the  first 
1.3  cm  downstream  from  the  apex.  This  modification  isolated  the  epoxy 
fiberglass  plate  form  the  discharge  and  prevented  further  ablation  from 
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Figure  21.  Comparison  of  Magnet  Impulse  and  Plasma 
Impulse  Measurements 


that  source.  At  the  same  time  the  sonic  slit  was  closed  and  gas  injection 
was  moved  to  a pair  of  1/8-inch  holes  in  the  upper  and  lower  electrodes. 

This  change  was  made  to  simplify  the  gas  supply  system.  Efficiencies 
greater  than  100%  were  again  observed  indicating  that  spurious  mass 
addition  must  still  be  occurring. 

To  obtain  a more  direct  location  of  the  plasma  current,  a miniature 
Rogowsky  probe  was  used. 

It  was  roughly  rectangular  In  shape  with  5 cm  long  sides  spaced 
0.5  cm  apart;  the  long  dimension  was  In  the  axial  direction.  An  axial 
scan  of  the  total  current  passing  through  the  probe  window  is  plotted  in 
Figure  22.  It  can  be  seen  that  no  current  was  running  forward  of  the  1 cm 
position,  l.e.,  the  measured  current  remained  constant  to  that  axial 
station.  It  then  decreased  rapidly,  falling  to  20%  at  2 cm.  This  indicated 
that  the  plasma  current  was  concentrated  in  a layer  a few  mm  wide  at  about 
1 cm  from  the  apex.  It  was  suspected  that  epoxy,  that  was  used  to  seal  the 
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slit- at  the  apex  (between  the  glass  side  plates),  was  supplying  mass  by 
ablation. 

This  region  was  covered  by  placing  another  glass  plate  against  the 
inside  surface  of  one  of  the  sides  and  moving  it  forward  to  close  off  the 
slit.  A plasma  current  scan  was  then  made  and  the  results  are  plotted  in 
Figure  22  for  comparison  with  the  uncovered  slit.  The  data  extrapolate  to 
about  1000  amperes  at  Z a 0.  The  measured  circuit  current  for  this  test 


Figure  22.  Axial  Variation  of  Current  in  Rogowsky  Probe 
(450  volts) 

was  3200  amperes.  Since  the  probe  area  covered  about  one  third  of  the 
current-carrying  cross  section,  the  probe  measurement  appeared  to  indicate 
that  the  current  density  was  roughly  constant. 

After  several  hundred  firings  of  the  glass  sidewall  thruster,  with  the 
modification  described.  It  was  found  that  the  glass  had  been  ablated  by 
discharges  originating  near  the  downstream  edges  of  the  large  brass  elec- 
trodes. This  showed  that  more  care  would  have  to  be  taken  to  reduce 
ablation,  and  work  on  the  nearly  completed  model  shown  in  Figure  23  was 
terminated.  This  thruster  was  built  with  a long,  narrow  plenum  to  permit 
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Figure  23.  Glass  Side-Wall  Axial  Injection  Thruster 
Designed  for  Reduced  Ablation 

axial  gas  injection  between  the  closely  spaced  magnet  pole  pieces.  But  it 
also  used  brass  upper  and  lower  electrode  surfaces. 

An  all -ceramic  thruster  channel  was  then  built  which  used  sonic  gas 
injection  through  copper  tubes  inserted  in  the  upper  and  lower  ceramic 
plates. 

The  initial  tests  with  this  model  showed  the  same  unrealistically 
large  Impulse  values  that  had  been  seen  with  the  glass  sidewall  thrusters. 
Because  of  the  clean  design,  this  was  not  thought  to  be  caused  by  ablation. 

A partial  explanation  appeared  to  be  the  interaction  between  the  current 
flowing  in  the  electrodes  and  current  leads  and  the  fringing  magnetic  fields 
above  and  below  the  magnet  pole  pieces.  This  possibility  was  examined 
experimentally  by  conducting  the  current  between  the  electrodes  in  a manner 
that  avoided  the  magnet  field.  This  was  done  by  attaching  copper  straps 
to  the  electrodes  Inside  the  channel  and  extending  them  downstream,  parallel 
to  the  thruster  axis  for  a distance  of  about  18  Inches  before  crossing  the 
electrode  gap. 


The  magnet  deflection  was  found  to  be  nearly  equal  to  that  observed 
in  normal  thruster  operation.  To  eliminate  this  effect,  the  current  was 
returned  through  a conductor  placed  between  the  pole  pieces  and  the  maqnet 
blocks  at  a location  such  that  no  magnet  swing  was  observed  when  the 
capacitors  were  discharged.  This  conductor  can  be  seen  in  Figure  11. 

The  thruster  was  then  operated  In  its  normal  mode  and  very  little 
magnet  swing  was  observed.  This  indicated  that  the  results  achieved  with 
the  glass  sidewall  model,  using  upper  and  lower  surface  gas  injection  may 
have  been  spurious.  Evidently  the  impulse  delivered  to  the  magnet  had  been 
due  to  the  interaction  of  the  lead  current  with  the  fringing  magnetic 
field  while  the  plasma  impulse  measured  by  the  ballistic  pendulum  inside  the 
vacuum  chamber  was  due  largely  to  the  ablated  mass.  Thus  these  impulse 
values  were  produced  by  different  mechanisms  (both  of  which  were  spurious) 
and  the  good  agreement  in  the  magnitude  of  the  effects  was  fortuitous. 

As  shown  in  Figures  24  and  25,  substantial  ablation  was  also 
encountered  with  this  model.  Figure  24  is  a close-up  view  of  one  of  the 
gas  injection/electrode  tubes  while  Figure  24  shows  the  build-up  of  ablated 
material  that  developed  on  all  four  of  the  channel  surfaces.  Because  of 
the  color  of  this  material  it  was  thought  to  be  copper  only. 


Figure  24.  Erosion  of  Copper  Gas  Inlet/Electrode  in 
Ceramic  Model 


Figure  25.  Ablated  Material  on  Inner  Surfaces  of  Ceramic  Model 
with  Copper  Gas  Inlet/Electrodes 

At  this  point  it  was  suspected  that  the  failure  of  the  channel  to 
produce  thrust  was  the  result  of  the  method  of  gas  injection.  The 
colliding  flow  between  the  hollow  electrodes  stagnated  at  hypersonic  speeds 
such  that  a pressure  of  several  torr  resulted.  This  pressure  was  too  high 
for  a low  resistance  discharge  so  that  the  current  would  tend  to  flow  well 
downstream  of  the  bias  magnet  field.  A self-field  acceleration  would 
occur  in  this  discharge,  but  no  impulse  would  be  delivered  to  the  magnet. 

To  overcome  this  difficulty,  and  hopefully  reduce  the  electrode  ablation, 
the  thruster  shown  in  Figure  11  was  built.  This  also  used  ceramic  walls 
for  all  surfaces  but  the  electrodes  were  solid  and  made  of  tungsten.  As 
with  the  first  model  tested,  gas  injection  was  through  a sonic  slit  about 
0.05  mm  wide. 

Again  large  magnet  impulse  measurements  were  made.  Indicating  more  than 
100%  efficiency  at  low,  but  reasonable,  argon  flow  rates.  After  about  100 
shots  the  ablation  pattern  shown  in  Figure  26  was  seen.  The  material  was 
black  and  was  thought  to  be  tungsten.  The  extensive  melting  of  the  anode 
shown  in  Figure  26  seemed  to  support  this  view.  One  additional  attempt  was 
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Figure  26.  Ablation  in  Ceramic  Thruster  with 
Tungsten  Electrodes 


Figure  27.  Tungsten  Anode  Ablation 
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made  to  stabilize  the  current  in  the  bias  field  region  by  placing  a 
ceramic  rod  between  the  electrodes  as  shown  in  Figure  28.  After  several 
discharges,  with  no  significant  change  in  performance,  the  ablation  shown 
in  Figure  29  was  found.  Not  shown  in  this  figure  is  a severe  erosion  of 
the  leading  edge  of  the  rod.  Examination  of  this  surface,  and  the  walls, 
with  a microscope  showed  that  much  of  the  ablated  material  was  alumina. 
The  ceramic  surfaces  resembled  a lava  flow. 

The  fact  that  large  magnet  impulse  values  were  measured  for  both  the 
open  channel  and  the  ceramic  rod  (flame  arrester)  configuration  indicates 
that  substantial  upstream  current  must  flow,  presumably  in  the  sidewall 
boundary  layers  and  carried  by  the  ablated  alumina.  The  ablation  pattern 
in  Figure  27  supports  this  interpretation. 
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Figure  29.  Ablation  of  Channel  and  Anode  After  Testing 
Ceramic  Rod  "Flame  Holder" 

With  this  result,  it  was  decided  that  the  high  power  densities  for 
this  mode  of  operation  preclude  the  development  of  a practical  thruster. 
One  further  test,  made  with  a mixture  of  nitrogen  and  hydrogen 
(N2  + 2 H2)  provided  dramatic  evidence  that  bias  field  thrust  augmentation 
is  not  feasible.  For  this  gas  mixture,  the  radiation  power  density  was 
evidently  too  low  to  produce  ablation  and  no  magnet  deflection  was 
observed. 


6.  CONCLUSIONS  • 


1.  Theoretical  calculations  of  efficiency  for  the  bias  field  thruster 
Indicate  large  performance  gains  relative  to  self-field  thrusters. 

2.  In  the  absence  of  ablation,  using  a mixture  of  nitrogen  and  hydrogen 
gases  (N,  + 2 Hz)  the  plasma  current  could  not  be  stabilized  in  the 
region  ot  the  bias  field  so  that  no  impulse  gain  due  to  the  bias  field 
was  obtained. 

3.  With  argon,  substantial  ablation  of  the  alumina  channel  and  tungsten 
electrodes  was  found  to  be  unavoidable. 

4.  It  is  concluded  that  a practical  bias  field  thruster  is  not  feasible. 
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